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Abstract —This work studies the performance of a cooperative 
network which consists of two channel-coded sources, multiple re¬ 
lays, and one destination. Due to spectral efficiency constraint, we 
assume only one time slot is dedicated for relaying. Conventional 
network coding based cooperation (NCC) selects the best relay 
which uses network coding to serve two sources simultaneously. 
The performance in terms of bit error rate (BER) of NCC, 
however, is not available in the literature. In this paper, we 
first derive the closed-form expression for the BER of NCC and 
analytically show that NCC always achieves diversity of order two 
regardless the number of available relays and the channel code. 
Secondly, motivated by a loss in diversity in NCC, we propose 
a novel relaying scheme based on partial relaying cooperation 
(PARC) in which two best relays are selected, each forwarding 
half of the codeword to help one source. Closed-form expression 
for BER and system diversity order of the proposed scheme are 
derived. Analytical results show that the diversity order of PARC 
is a function of the operating signal-to-noise ratio (SNR) and 
the minimum distance of the channel code. More importantly, 
full diversity order in PARC can be achieved for practically 
operating finite SNRs with the proper channel code. Elnally, 
intensive simulations present a huge SNR gain of PARC over 
NCC and reference schemes without relay selection. 

Index Terms —Cooperative diversity, relay selection, partial 
relaying, channel coding. 


I. Introduction 

In wireless networks, idle nodes have a potential to par¬ 
ticipate in transmission of other nodes to form cooperative 
communication. Cooperation among nodes has been shown as 
an effective technique to widen the coverage and to improve 
the performance of wireless networks in both terms of Signal- 
to-Noise Ratio (SNR) and diversity gain Q. In the most basic 
cooperation form of single-source single-destination network, 
a relay estimates the source signal and then forwards it to 
the destination. It is shown that relay networks achieve a 
performance gain when compared with the non-cooperative 
counterpart Q. In order to achieve this gain, however, addi¬ 
tional orthogonal channel is usually required, which results in 
loses in spectral efficiency, especially when more than one 
relay is employed. Fortunately, such loss in multiple-relay 
networks can be effectively reduced by using opportunistic 
Relay Selection (RS), in which only the best relay is selected 
for cooperation Q. It is shown that RS can achieve full 
diversity order for single-source multiple-relay networks and 
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outperform other relaying schemes in terms of SNR gain and 
effective capacity 0 - 

Network Coding (NC) has gained tremendous attention 
because of its potential improvement in diversity gain and 
throughput over classical routing techniques |[^. In NC, an 
intermediate node combines multiple input packets into a 
linear combination which is then forwarded. Recently, there 
have been much studies on combining NC together with RS 
to further improve the spectral efficiency, mostly focusing 
on the Two-Way Relay Channel (TWRC). The authors in 
Q propose a joint design of NC with RS for Decode-and- 
Forward (DF) TWRC based on the max-min criterion in order 
to maximize the worst relay channel. In g, a SNR-based 
suboptimal relay ordering is proposed for two way Amplify- 
and-Forward (AF) relay networks. A similar method is studied 
in Q to derive the system Outage Probability (OP), BER, and 
diversity order. Compared with research on RS in TWRC, 
which shows full diversity is achieved and is frequently 
available in the literature, research on RS in unidirectional 
relay networks is still limited. The study of NC with RS 
in unidirectional networks is first considered in 0 - In this 
work the authors study the Diversity Multiplexing Tradeoff 
(DMT) and show that full diversity is achieved. However, 
the analysis in 0 is relied on an unrealistic assumption that 
unintended packets are available at all destinations, which 
simplifies the unidirectional networks to TWRC. A generalized 
DMT analysis is presented in flO) . Likewise, the authors in 
1101 also assume an optimistic assumption that the selected 
channels are independent, whereas these channels belong to 
an ordered SNR sequence and hence are highly correlated 
Q. The analysis of the counterpart AF relaying in inter-user 
interference channels are studied in UTj-jn). We note that the 
diversity order in the above-mentioned works is studied via the 
limit of the upper bound of either OP or BER. Eurthermore, 
most of these papers did not consider channel coding, which 
might be in contrast to practical scenarios in which nodes are 
usually protected by some forward error correction codes. 

In this paper, we investigate the performance of coopera¬ 
tive networks under practical conditions, i.e., the transmitted 
signals are protected by Convolutional Codes (CC). In the 
considered system, two sources communicate with a com¬ 
mon destination with the aid of multiple available relays. 
This scenario can find applications in uplinks cellular mobile 
systems where two mobile users try to send data to the base 
station and some surrounding friendly, idle users can act as the 
relays. Due to the spectral efficiency constraint and processing 
delay limit, it is assumed that only one timeslot is dedicated 
for cooperation. The best RS is employed Q to effectively 
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exploit the spatial diversity. For low-complexity functions at 
the relays, Demodulate-and-Forward (DMF) relaying protocol 
0 is used. At the destination. Cooperative Maximal Ratio 
Combining (C-MRC) detector GD is used priori to channel 
decoding to avoid error propagation. It is worth to note that 
C-MRC is a suboptimal detector and provides full diversity 
gain and a performance close to Maximum Likelihood (ML) 
receiver ID- 

The contributions of the paper are as follows: 

• Firstly, we investigate the performance of Network Cod¬ 
ing based Cooperation (NCC) in which one selected relay 
helps two sources by applying network coding on the 
estimated codewords. Closed-form expression of the BER 
is derived for all sources. From the analyzed BER, we 
analytically show that NCC always achieves a diversity 
of order two regardless the channel code and the total 
number of relays. This result coincides with the diversity 
order derived from OP analysis 0, 0. 

• Secondly, we propose a novel relaying scheme named 
Partial Relaying based Cooperation (PARC). The key 
difference between PARC and NCC is that in the former 
two relays are selected for cooperation, each helping one 
source. Due to the spectral efficiency constraint, each 
selected relay in PARC forwards half of the estimated 
codeword to the destination. The cooperation based on 
partial relaying has been studied by some authors in 
0, 0. Compared with these works, our proposed 
scheme has two main differences: i) we investigate the 
system via BER analysis, whereas these papers study the 
system through OP, which is fundamentally different from 
our method; and ii) the proposed PARC employs RS to 
improve the spectral efficiency, while these papers do not. 

• Thirdly, insightful theoretical analysis is provided for 
PARC in finite SNR regime. Particularly, closed-form 
expression for BER and the diversity order are derived, 
which reveals that the instantaneous diversity order of 
PARC is a function of the operating SNR and the min¬ 
imum distance of the channel cod^ More importantly, 
PARC can achieve full diversity order with suitable chan¬ 
nel codes in low and mediate SNR regime, the operating 
SNR region in practical systems. Intensive simulation 
results show a large SNR gain of PARC over NCC and 
other reference schemes. 


The rest of the paper is organized as follows. Section III 
describes the system model of PARC and NCC. Section [III 
provides details for the relay selection process. Section ^ 
analyzes BER and diversity order of PARC. The performance 
analysis of NCC is analyzed in Section jV] Section jVT] shows 
numerical results. Finally, conclusions and discussions are 


given in Section VII 


IT System Model 

The system under consideration consists of two sources de¬ 
noted by Si and S 2 , Nr relays denoted by with 1 < * < Nr, 

*To investigate the system in finite SNR regime, instantaneous diversity 
order is defined as the generalized definition of classical diversity order at 
any SNR value. More details are presented in |IV-B| 


and one destination denoted by D. All nodes are equipped with 
a half-duplex single antenna. The system is assumed to operate 
on orthogonal channels with perfect time synchronization. As 
a result, a cooperation period is divided into two phases: broad¬ 
cast phase and relaying phase. Due to the spectral efficiency 
constraint and processing time limit, we assume that only one 
time slot is dedicated for relaying. All the channels are subject 
to block Rayleigh fading plus Additive White Gaussian Noise 
(AWGN). In order to minimize the computation at the relays, 
Demodulate-and-Forward (DMF) relaying protocol is used. To 
effectively exploit spatial diversity and to achieve high spectral 
efficiency, relay selection is used 0. The relay selection 
process is performed at the beginning of every cooperation 


period and will be described in details in Section III 


A. Partial Relaying based Cooperation (PARC) 

Motivated by our previous work which showed full diversity 
gain is achieved for three-node relay networks in low and 
medium SNRs even when the relay only forward parts of 
the codeword 0, we propose to combine relay selection 
and partial relaying. In the proposed PARC, two relays are 
selected in which each is the best relay for one source. Since 
two relays are active in the relaying phase, each relay only 
occupies half of relaying time slot, as shown in Fig. [T^. 
Consequently, the selected relay can only forward half of the 
estimated codeword to the destination. Although forwarding 
half of the estimated codeword, full diversity gain in finite 
SNR^is expected when suitable channel code is employed, 
as shown later on in Section ITVl 

First, source Si,i = 1,2, encodes a AT-length data message 
Uj into a codeword which contains N coded symbols by 
a convolutional code g with code rate K/N. The codeword 
Ci is then modulated into a signal x^. Next, the signal is 
broadcasted to the relays and the destination. Denote , Rs^ 
as the selected relays for Si and S 2 , respectively. The received 
signal at the destination and the selected relays at the end of 
first phase are given as follows: 

f ySiRs, = s/PSiRs,hSiRs,^i + ns,fis. = 1,2, 

1 ySiD = s/PSiOhSiD^i + TtlSiD,i = 1,2, 

where Pxy with X S {S'i,S' 2 },V € {i?s^, D} is the 
average received power at node Y from node X, including 
the path loss; hxY is the channel fading coefficient between 
X and Y, which is a complex Gaussian random variable with 
zero mean and unit variance, i.e., E||/ixyp} = 1, and is 
mutually independent among X ^ Y channels; n() is a noise 
vector whose components are Gaussian random variables with 
mean zero and variance a^. 

At the end of the first phase, the selected relay estimates the 
source coded symbols and forwards them to the destination. 
In the proposed scheme, the selected relay Rs.,i = 1,2, 
uses half of the relaying time slot to forward half of the 
codeword to the destination. More specifically, the selected 
relay Rs^ first estimates L = [1V/2J source coded symbols to 
form an estimated punctured codeword 

^It is also called instantaneous diversity order, which is measured as the 
slope of the BER curve plotted in log-log scale as a function of SNR. 
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Fig. 1. Time allocation for the Partial Relaying based Cooperation (a) and Network Coding based Cooperation (b). In PARC, two relays are selected out 
of Nr total number of available relays, each forwarding half of source codeword. On the other hand, only one relay is selected in NCC. The selected relay 
RiVC ia NCC forwards the whole network-coded codeword to help two sources simultaneously. 


where [aj denotes the largest integer less than a, and 0 = 
{fci, k 2 , ■ ■ ■ ,kL} being the set of indexes of the symbols 
which are helped by the relay R 5 .. The index set 0 are 
determined randomljjj The source coded symbols at the relay 
are estimated by the ML detector as follows: 

CRs^,i = arg min {\ys,Rs,M “ 

Vki G 0 , where Xi^ki being the corresponding modulated 
symbol of Ci^ki- Next, modulates into the modulated 
signal x/jg and then forwards it along with the index set 
0 to the destination. The cost for conveying the index set 
is negligible since it can send, e.g., the seed of the random 
interleaver, to the destination. 

The received signal at the destination transmitted from the 
relay is given as: 


YRs^D 



Rs-D^Rs^ + nRs-D,i = 


1 , 2 , 


( 2 ) 


where u is the channel coefficient from Rg. —> D, and 
^Rs D is a noise vector whose components are Gaussian 
random variable with zero mean and variance cr^. 

After receiving two signals from the source and the relay, 
the destination starts the decoding process with two consec¬ 
utive steps: demodulating and decoding. Assuming that the 
CSl of all channels, i.e.. Si —D,Si —>^ Rs- and Rs. —>^ D 
channels, are available at the destination, the destination first 
applies the C-MRC detector HD to demodulate the coded bits 
for source Si, i = 1 , 2 , as follows: 

Ci k = arg min M{xi^k), ^ < k < N, 

Ci. fcG{0,l} 

where the detection metric M(xi^k) = \ySiD,k — 
\/PSiDhsiDXi^k\^ if fc ^ 0 ; otherwise 


M{Xi^k) = ySiD,k — \J Ps^ohSiDXi^k 


(3) 


+ A 




yRsiD,k — PRs^ohRs^DXRs.^k 


In ([^, Ai^g is the parameter of the C-MRC detector which 

’ A min(7SiHs .7Hs d) 

IS computed as Aijg. = - —-—, where ^xy = 

being the instantaneous SNR of the channel 

X ^ y. 

The C-MRC detector then computes log-likelihood ratio 
values of the coded bits and sends them to the channel decoder. 


^Other selection of ©, e.g., optimal index set, can be employed, but are 
beyond the scope of this paper. 


Finally, the channel decoder applies the BCJR algorithm pT| 
to decode the transmitted data. 

Remark 1: In our protocol, the selected relay always for¬ 
wards the estimated symbols to the destination. Fortunately, 
possible decoding error in c/jg hence error propagation, is 
effectively mitigated by Xr^. in C-MRC. For example, if the 
source-relay channel is too noisy, i.e., 'jSiRs- too small, it 
is highly probable that R 5 . decodes with errors. In this case, 
however, Xs^ is small and the contribution of the relayed signal 
is negligible. 


B. Network Coding based Cooperation (NCC) 

In NCC, the relays use network coding to help both sources 
simultaneously to improve the spectral efficiency. One cooper¬ 
ation in NCC is also divided into two phases: broadcast phase 
and relaying phase. The broadcast phase is similar to that in 
PARC, whereas in the relaying phase, only one best relay is 
active. Time allocation of NCC is depicted in Fig. [J. Unlike 
PARC, the selected relay in NCC forwards the whole network- 
coded codeword to the destination. Without loss of generality, 
denote by Ratc the selected relay in NCC. The received signal 
at Ratc is given as follows: 

YSiR^c sj PSiRwc kSi Rnc “f ^SiRNc U (^) 

At the end of the first phase, Knc decodes the estimate Xifj 
of Xi,i = 1, 2, using the ML detector as follows: 

CiR,k = arg min “ V^SiRNc^SiRivc^iM"^}^ 

Ci, feG{0,l} 

* e {1,2}, 1 < fc < TV, 

where Xi^k being the corresponding modulated symbol of Ci^k- 
Then Knc performs network encoding to get cxc, where 
CNC,k = Cift k ® C 2 R,k, Vfc, and 0 denotes the binary XOR 
operation. 

The received signal at the destination from the selected relay 
is given by: 

YRncD = s/ PrncD^RncD^nc + (5) 

where ixc is the modulated signal of cxc- After two phases, 
the destination receives three channel observations from two 
sources and the selected relay. To decode the source data, the 
destination applies joint network/channel decoding algorithm 
to a ’’compound code” G pO) which sees the relayed signal 
as additional parity bits (redundancy). The compound code G 
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is formed from the individual code g as follows: 


G = 


g 0 g 
0 g g 


( 6 ) 


where 0 is a zero matrix withs same size as g. For full details 
of joint decoding at the destination, we refer the readers to 



III. Relay Selection FOR PARC AND NCC 

In this section, we describe in details the relay selection 
processes for PARC and NCC and provide essential charac¬ 
teristics of the selected relay channels. 

A. Relay Selection in PARC 

The relay selection is based on the suboptimal max-min cri¬ 
terion that maximizes the worst end-to-end SNR and reduces 
computational complexity Q. The relay selection process in 
PARC is executed for each source separately and can be 
done in distributed manner similar to |(3). After the channel 
estimation, the relays set a timer that is inversely proportional 
to their channel gain. The first relay whose timer is zero will 
send a pulse to the destination. Upon receiving the pulse, the 
destination declares the chosen relay. Because the selection 
procedure for two sources are similar, we avoid source sub¬ 
script in this subsection for notation brevity. Particularly, the 
source is denoted by S and the selected relay is denoted by 
Rs. Since the relayed symbols received at the destination 
using DMF protocol can be well described by equivalent 
channel we model a two-hop source-relay-destination link 
by an equivalent single-hop channel 7 ^ = min { 75 / 2 ^, 7 ^^^ £>}, 
1 ^ j ^ A^r- In Rayleigh fading channel, both ysRj 

exponential random variables with mean 75 ^ and 
7 ^ £), respectively. Using the property of the Min function, it 
is straightforward to show that 7 _, is also an exponent random 
variable with mean 7 ^, which is computed as: 

1 _ 1 1 

Ij ISRj ^RjD 

To minimize possible errors of the relayed symbols, the 
relay that has the biggest equivalent channel is selected for 
cooperation: 

Rs = argmax 7 j, 1 < j < Nr- 
The equivalent channel of the selected relay, ysei, is given by: 

ISei = max{7i,... ,7ArU- 

By using the Max function | [22| , the Probability Density 
Function (PDF) of ^Sei is given in a shorten form as follows: 

Nr- 

ni = l,...,nj = 



fls.i ( 7 ) = 51 
i=i 


1 


1 ^Sel,j 


exp (- 


‘y S el, j 


The moment generating function (MGF) of ysei is given 
by: 

Nr 

i-v’-' E 

ni=l,...,nj = 

Remark 2: In the proposed PARC, the relay selection is 
performed for each source separately. Also the decoding at 
the destination is executed separately for each source. 


ll + 7SeijS 


(7) 




Nr 

..w = E( 


i=i 


B. Relay Selection in NCC 

The relay selection process is performed in NCC based 
on a criterion that minimize possible error of network-coded 
symbols. Because an error of the network-coded signal can 
result from either source-relay links or relay-destination link, 
the network-coded symbols can be seen as if it has been trans¬ 
mitted via an equivalent channel which yields the same error 
probability |T§. Using the equivalent error probability for 
network-coded symbols, the two-hop source-relay-destination 
channel corresponding to the relay R^ can be tightly modeled 
as follows m- 

leq,j = min{7Sifl,, -is-,R,, 1 r,d}- 

Because ysiRjTlS^Rr ^re exponential random 

variables, it is straightforward to show that yeq,j is also an 
exponential random variable whose mean j is given by 


leq,j ISiRj TSsflj 7fljD 

The best relay in NCC, denoted by Knc- is selected by the 
max-min criterion as: 


Rac = argmax7eq7, 1 < j < A^r- 
Rj 

The equivalent network-coded channel of the selected relay is 
chosen as follows: 


JNC = max{ 7 eq,i, . . . ,7eg.AU- 

Because the 'yeq,j are mutually independent, the Cumulative 
Density Function (CDF) of yjqc is computed as: ( 7 ) = 

n^i j ( 7 )- Taking the derivative of Fiyjvc(7) obtain 
the PDF of jNc expressed in the simplified form as follows: 


fjNC ( 7 ) = 


1 


Nr Nr 

j=l ni = l,...,nj=l 7aCj' 


where 


rij , 


1 1 

-f 


exp - 


1 


^ 7AC,i ^ ^ 


ISiRk "IS^Rk ^RkD 


The MGF of is calculated as follows: 

Nr Nr 


1 


i=i 


ni = l,...,71^=1 




(8) 


where 
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IV. Performance Analysis for Partial Relaying 
BASED Cooperation 

In this section, we analyze BER and diversity order of PARC 
using the equivalent channel model. Because the decoding pro¬ 
cedure for two sources is similar, the analysis for two sources 
is analogy. To simplify notation in the analysis, we drop source 
subscript. After two phases, the destination receives two signal 
from S and R 5 , the selected relay. The combined signal at the 
output of C-MRC detector can be classified into two groups: 
the first group consists of symbols which are helped by the 
selected relay, and the second group includes the rest symbols 
which are not relayed. In other words, the received signal at the 
destination can be seen as an output of block fading channel 
with 2 blocks: one block consisting of N — L symbols only 
sees the channel ■jsd, and the other one which contains L 
symbols sees both channel jsd and channel jsei- 


A. Bit Error Rate Analysis 

Let Pu{d) be the Unconditioned Pair-wise Error Probability 
(UPEPj^ which is probability that the destination decodes for 
a codeword with the Hamming distance d (number of non-zero 
coded bits in c^) when the all-zero codeword was transmitted. 
The BER of PARC is upper-bounded as follows | |2^ : 

N 

Pe < E w{d)Pu{d), (9) 

d=f 

where / is the minimum distance of the channel code g, 
and w{d) is input weights which is number of non-zero 
information bits in U; and is computed directly from structure 
of the code. The UPEP Pu{d) is the expectation over the 
channel fading coefficients of the Conditioned Pair-wise Error 
Probability (CPEP) Pc(d): Pu(d) = E{Pc(d)}. The CPEP 
Pc((i) obviously depends on the fading channels and how d 
non-zero coded bits are distributed on the two blocks (jsd 
and JSD + ISei)- Denote = {di,d2}, di + d 2 = d, as 
the weight pattern that presents how d weights are distributed 
on the two blocks. Because d non-zero coded bits uniformly 
locate in the two blocks, the CPEP can be further analyzed as 
follows: 


Pc(d) = 5 ]Pc(d|Dd)p(Dd), ( 10 ) 

Dd 


where p(Dd) is the probability of pattern D^, which is 
computed by combinatoric computation as 


p(Dd) 


pN ■ 


where 

Substituting (|T0|) into Pu{d) we obtain: 


P,(d) =^E{Pc(d|Dd)}p(Drf). (11) 

P„(d|Drf) 


Given the pattern = {di, d 2 }, there are di non-zero coded 
bits undergoing through block js d and d 2 non-zero coded bits 


“^The unconditioned pair-wise error probability does not depend on the 
fading channels. 


undergoing through block jsd + 7Se/- As a result, the CPEP 
Pc((i|D(i) is calculated, by using the same techniques in | [24) , 
as follows: 

Pc(ti|Dd) = Q (v^) , (12) 

where j^ = dijsD +^ 2(750 + ISei) = djsD + d 2 jsei and 
Q{x) = ^"^dt denotes the Q-function. 

Taking into account the independence between jsd and 
jSeU we obtain the UPEP Pu(d|D£;) given in Theorem 
below. 

Theorem 1: Given the weight pattern = {di,d 2 },rf = 
di -f d 2 , the UPEP Pu ((i|Dd) of PARC is given as follows: 


' 1 (> - . if <i 2=0 


Pu{d\T>d)={ E((-1V ^ E d:i{djsD,d2JSel.j)), 
if fi 2 > 0 


where 


Proof: See Appendix [A| ■ 

Substituting UPEP in Theorem into o and (|^ we obtain 
the upper bound for the BER. Note that even though d in (|^ 
can be as large as the codeword’s length, i.e., N, the BER 
usually depends on few first values in fading channels. To 
give insightful understanding of PARC, we analyze the system 
diversity order. 


B. Finite-SNR Diversity Analysis 

The classical definition of diversity order is defined as the 
negative exponent of the average BER as a function of SNR 
in log-log scale when the SNR tends to infinity | |25| , which 
visually represents the slope of BER curve in infinity SNR 
domain. In this paper, since we are interested in finite SNR 
regime, we define the diversity order at a certain SNR 7 as 
follows: 


C(7) = - 


lim 

A-)-l 


log[P2(A7)]-log[Pe(7)] 

log(A 7 ) - log( 7 ) 


(13) 


which obviously matches the classical definition of diversity 
when the SNR tends to infinity. Because the diversity order 
depends on the average SNR, we refer to (^( 7 ) as instantaneous 
diversity order. The key idea behind the definition is that it 
allows to study the behavior of the system at any SNR values. 

We might write x = SAR”** if x has asymptotic diversity 
order rj (classical definition of diversity order). Erom Q we 
know that the diversity order of PARC is determined by Pu{d) 
because the input weight w{d) of the channel code is constant. 
We first compute diversity order of the UPEP for a given 
weight pattern as in Theorem 

Theorem 2: Given the weight pattern = {di,d 2 } with 
d = di+d 2 , the UPEP Pu (d|Dd) in PARC is given as follows: 


Pu id\T>d) A 


SNR-\ if(i 2=0 

SNR-^^^+P, ifd 2>0 
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Proof: See Appendix [B] ■ 

Theorem states that Pu(d|D(i) can have either diversity 
order one or diversity order AV + 1. From o we conclude 
that Pu{d) is a combination of one factor with diversity of 
order one and one factor with diversity of order + 1- The 
contribution of the factor with diversity order equal to one is 
computed from ([T4ll as follows: 


^NI2 

p{D,^{d,0}) = ^ 


N-2k 
-l-J- 2N-2k' 

k=0 


In practical systems, the codeword length N is usually much 
larger than d, then p{Di) can be well-approximated as 

Substituting the result in Theorem into we can write: 

Pe = A:2-^SNR-^-f SNR-^’--^ (15) 


where K is the normalized constant that depends on the 
channel code and network topology. 

From •El and EH* we conclude that the instantaneous 
diversity order of PARC consists of one factor which achieves 
full diversity order and one factor which achieves a diversity 
of order one. The impact of the diversity order one factor is 
inversely proportional to the channel code’ strength, i.e., its 
minimum distance. For a strong code with large minimum 
distance /, the impact of the diversity one factor is negligible 
compared with the diversity order A^^ + 1 factor, resulting the 
system can achieve full diversity Af^ + 1 in the finite SNR 
regime. This result is important because the operating SNR 
regime in practical systems is usually finite. 


V. Performance Analysis eor Network Coding 
BASED Cooperation 

In this section, we analyze the BER and diversity order 
of NCC. Using the equivalent channel, the two-hop network- 
coded signal can be modeled as if it was conveyed by a single 
channel whose instantaneous SNR is 7 nc ED- 


0 ). To derive ( [T6] l, it requires to know the minimum distance 
F of the compound code, the input weight Wi{d) and how 
d bits in the compound codeword X are distributed among 
three channels Si —> D, S 2 —> D, and Ratc —> D. Denote 
Wd = {di,d 2 ,dii} as the weight pattern that specifies how 
d weights are distributed among these channels, where di is 
the output weight of the individual codeword transmitted via 
the channel ^ D or channel Ratc —^ D. By definition, 
d — di + d 2 + da- The input weight and the pattern can 
be computed via heuristic searching of the trellis of G. The 
following result is important for further analysis. 


Lemma 1: The minimum distance F of the compound code 
G is equal to twice the minimum distance of the single code g, 
F = 2/, and the weight pattern Wp has one of the following 
values {/,/,0}, {/,0,/}, {0,/,/}. 


Lemma 2: For any pattern of the 

compound codeword X with output weight d > F, there are 
at least two non-zero elements in W^. 

The proof of Lemma [Tl and Lemma are given in | [20) . 
Lemma and Lemma^ provide an important information 
about the output weights of the compound code: d weights of 
the compound code always experience at least two independent 
channels. Lurthermore, the number of patterns is finite and 
strictly defined by G. 


Using Lemma [T] and 1^ we can reformulate ([T^ as follows: 

_ , +00 ^ 

d=F Wd 

where Pu((i|Wrf) is the UPEP depending on the pattern 
and is the expectation of the CPEP over the fading channels: 

Pu(li|Wd)=E{Pe(d|Wd)}. 


A. Derivation of BER 

Recalling that in NCC, the destination applies the BCJR 
algorithm on the compound code G, which is described in 
Section [ri-B| The compound code G has compound input X = 
[xi,X 2 ,XArc'] and channel output Y = [ysiU, yssD, yi 2 „cD]- 
Note that the output of G undergoes some block fading 
channels with three blocks 'ysiD,'ys 2 D, and jnc, and it 
decodes the data messages instantaneously. Consider G as a 
regular channel code, the BER of source Si is calculated as: 

_ , +00 

Pe^ = - ^ Wi{d)Puid), i = 1,2, (16) 

^ d^F 

where F is the minimum distance of the compound code G, 
Wi{d) denotes input weights corresponding to source Si in the 
compound codeword, and Pu(<i) is the UPEP of receiving a 
super codeword with the output weight d, assuming that the 
all-zero compound codeword has been transmitted (ci = C 2 = 


It is assumed that the erroneous detected symbol could only be 
one of the nearest neighbor symbols. Using the Gray mapping, 
each closest symbol error only causes one coded bit error. 
Therefore, the CPEP Pc(d|D(i) is approximated as pO) : 

Pc (d|Wrf) = Q (v/2^) , (18) 

where 72 ^^ = dijSiD + ^ 2752 ^ + is defined as 

the total SNR at the destination in NCC. Because the three 
channels in mutually independent, the MGE of the 

total SNR can be computed as follows: 

Applying the MGE method we can derive the UPEP 
Pu (d|Wj;) in NCC as in Theorem]^ 

Theorem 3: Given the weight pattern = {di,d 2 ,dfl}, 
d = di + d 2 + dn, the UPEP Pu ((i|Wd) of the compound 
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code in NCC has a form given by; 


' Ii (di7s^£,,d27s2-D) > ifdfl = 0 

ifdi=0 
j=i 


Pu (d|Wd) - .( ^ 


if d2 = 0 
if did2dji y 0 
where Xi (a, b) has been defined in Theorem 


f=i 

f=l 


2:1 (d27s^i),di?7jvcj). 


ni=l,...,nj=l 
Tlx^" '^'^j 


fin _ 




ni—l,...,nj—l 


fia — 


y]] I 2 (dl7SiD> d2-fs^D,dRjNC,3) > 


ni=l,...,nj=l 

Til 


and 


2:2 


fa, 6 ,c) = iff - 7 -- 

V / 2 V (a — 5)(a — c) Y o + 1 


62 

1 b C2 

/ ^ ^ 

(6- a){b-c) \ 

6-1-1 (c — a)(c — 6) \ 

' c-f 1 / 


Proof: See Appendix ■ 

It is worthnoting in ( [T7| l that the BER of each source in 
NCC is a sum of terms given in Theorem]^ weighted by their 
corresponding input weights (W^). In NCC, the weight 
pattern only holds a few values and the input weights are 
computed from the extended distance spectrum. Table [I] gives 
an example of the distance spectrum of the compound code. 


TABLE I 

Input weight and output weight distribution at d = F = 24 oe 
COMPOUND CODE G IN 0. g = [23, 35, 37] 


Wi 

W2 

di 

C?2 

dR 

0 

12 

0 

12 

12 

12 

0 

12 

0 

12 

12 

12 

12 

12 

0 


B. Diversity analysis 

Since the BER in NCC is a linearly proportional to the 
UPEP Pu (c?|Wd) via the input weights, the diversity order 
of NCC is equal to diversity order of the UPEP. Employing 
asymptotic equivalent notation as in the previous section, the 
diversity order of the UPEP is given as the below theorem. 

Theorem 4: Given the weight pattern = {di,(i 2 , 
with d = di+ d 2 + da, the UPEP Pu (fi| W^) in NCC has an 


asymptotic equivalent form as follows: 

r SNR’Y ifdit^O 

Pu {d\Wd) = I SNR-(^'-+i\ if di = 0 or d 2 = 0 . 

[ SNR-(^'-+^\ ifdid2dRy0 

The proof of Theorem ^ is given in Appendix 

It is shown from ([Hi and Theoremthat the BER in NCC 
is a combination of three factors whose diversity orders are 
respectively 2, + 1 and Nr + 2. As the contribution of 

these factors are comparable and equal input weights of the 
compound code (shown in Tablej^as an example), the diversity 
order of NCC is dominated by the diversity order 2 factor. 
Consequently, NCC achieves diversity order 2 regardless the 
channel code and the total number of available relays. 

VI. Numerical Results 

This section shows simulation results to confirm the ef¬ 
fectiveness of the proposed system described in Section |I^ 
All channels are subject to quasi-static block Rayleigh fading 
plus AWGN. Because we focus on the diversity order, and the 
modulation order does not change the system diversity order, 
BPSK modulation and binary network coding are carried out 
in simulations. The data packet length is equal to 1024bits. 
We consider symmetric network, i.e., fsiRj — ISR^lRjD = 
iRD^TSiD — The relays locate at the middle of 

the sources and the destination and the pathloss exponent 
is 3.5, resulting in that the average SNR in source-relay 
channels and relay-destination channels are 10.5dB better than 
source-destination channels. Note that our analysis validates 
for arbitrary locations of the relays. The channel code is 
chosen as the one that optimizes both the minimum distance 
and distance spectrum in block Rayleigh fading channels | [26) . 
Different channel codes g are compared. 

We also present the performance of two reference schemes. 
The first reference scheme (named Reference 1 in the figures) 
is based on fractional repetition coding cooperation (H, fig. 
The second reference scheme employs factional repetition 
coding together with network coding (named Reference 2 in 
the figures). All relays are active and share the relaying phase 
in two reference schemes. In Reference 1, since the relays 
help the sources separately, each relay forwards l/{ 2 Nr) of 
the estimated codeword. In Reference 2, all relays use NC to 
help the sources and each relay randomly forwards 1/Nr of 
the network-coded codeword. We note that no relay selection 
is used in the reference schemes. 

Eigure [^compares the performance of PARC and referenced 
schemes for the channel code [133 165 171] with code rate 1/3 
and the minimum distance / = 15. The total number of relays 
Nr equal to 2 and 3 are plotted. The observed performance 
region satisfies BER > 10“® because this is the target BER 
for most practical applications. It is shown in the figure that 
the proposed PARC achieves an instantaneous diversity order 
of 3 and 4 (full diversity order) in the observing SNR range 
when the total number of relays is 2 and 3, respectively. 
Such expected result can be explained from Theorem that 
in this case, the impact of the diversity one factor equals 
p(Di) = (1/2)-^ ~ 3.10“®, which is negligible. Therefore, 
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(a) Nr = 2 (b) Nr = 3 

Fig. 2. Performance comparison between PARC and NCC when the CC [133 165 171] with the minimum distance / = 15 and the rate 1/3 is used. 




(a) Nr = 2 (b) Nr = 3 

Fig. 3. Performance comparison between PARC and NCC when the CC [25 33 37] with the minimum distance / = 12 and the rate 1/3 is used. 


the diversity order of PARC is determined by the full diversity 
factor in the observing SNR region. In contrast, NCC always 
achieves a diversity of order 2, which is a consequence of 
Theorem A huge SNR gain is achieved by PARC. In 
particular, PARC outperforms all other schemes about 5dB 
for Nr = 2 and 9dB for = 3 at BER = 10"^. Another 
observation is that Reference 1 also surpasses NCC because 
the relayed symbols in Reference 1 see more spatial diversity 
gain than that in NCC. When SNR tends to infinity, NCC may 
outperform PARC because the diversity order of PARC will 
collapse to one while NCC still has diversity order equal to 
2. From the practical system point of view, this crossing-point 
might not weaken the advantage of PARC over NCC since 
practical systems usually operate at finite SNR regime. 

Figure shows the performance comparisons when the 
channel code [25 33 37] with rate 1/3 is used. The minimum 
distance of this code is equal to 12. It is not surprised that 


NCC always achieves diversity order 2 for both Nr = 2 or 
Nr = 3 and the performance of NCC in both cases is similar. 
It is observed that PARC only achieves full diversity order in 
low SNRs. More specifically, PARC achieves diversity order 3 
in the SNR range until lOdB when Nr = 2 and diversity order 
4 until SNR = 5dB when Nr = 3. When SNR increases, a 
degradation in instantaneous diversity order is observed, which 
is predicted by our analysis (for this code, the contribution of 
diversity order one factor approximately is 2~f ~ 2.4e — 4). 
However, a similar SNR gain as for strong code [133 165 
171] is achieved by PARC at BER of 10“^, which is about 
5dB for Nr = 2 and 7dB for Nr = 3. A sound interesting 
observation is that the performance Reference 1 for Nr = 3 
is worse than that for Nr = 2. This is because in the later, the 
relay forwards less symbols in W,, = 3 than in = 2 case, 
and the channel code is not strong enough to compensate for 
less relayed symbols in A^^ = 3 case |[24). 
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Fig. 4. Performance comparison between PARC and NCC for different 
minimum distances, which corresponds to code’s correction capacity. Three 
codes with rate 1/3 are compared: CC [5 7 5] with f = 7, CC [25 33 37] 
with / = 12, and CC [133 165 171] with / = 15. 


Figure compares the performance between PARC and 
NCC for different channel codes and Nr = 3. Three codes 
with different error correction capabilities are presented: the 
weak code [5 7 5] with small minimum distance / = 7, 
the moderate code [25 33 37] with / = 12, and the strong 
code [133 165 171] with / = 15. Full diversity order is 
observed in low SNRs for all codes. When SNR increases only 
the strong code achieve full diversity order. The weak code 
starts losing diversity order earliest at SNR of 5dB, while the 
moderate code’s diversity degrades at lOdB. Compared with 
NCC, however, PARC significantly outperforms for all codes 
in the observed SNRs. 

In conclusion, the most effectiveness of the proposed PARC 
is that it can achieve full (instantaneous) diversity order in the 
low and medium SNR regime, which in turn results in a large 
SNR gain in finite SNR region. This is crucial for practical 
systems because their operating SNRs is usually finite. 

VII. Conclusions and Discussions 

We have proposed a novel cooperative scheme for a two- 
source multiple-relay networks that combines relay selection 
and partial relaying to effectively exploit the spatial diversity. 
In the proposed scheme, the selected relay retransmits half 
of the estimated codeword in order to satisfy the spectral 
efficiency constraint. We have analytically shown that the 
proposed scheme can gain full diversity order in finite SNR 
regime (instantaneous diversity) when a suitable channel code 
is used. It has been shown that the instantaneous diversity 
order is a function of the minimum distance of the code and 
the operating SNR. Numerical results show a significant SNR 
improvement the proposed scheme compared with reference 
schemes. 

The proposed partial relaying cooperation can easily be 
extended to general networks that consist of multiple sources 
and a given finite time slots for relaying. In this case, the 
selected relays might forward a number of symbols different 
from half of codeword. A promising application of PARC is 


to design for the sources which have different error correction 
capacities to achieve a given target BER. As such, how many 
relayed symbols for a source should be carefully chosen 
depending on the strength of its channel code. 


Appendix A 
Proof of Theorem[T] 

Because the relayed symbols are randomly distributed on 
the codeword, the weight d 2 on the relayed block can have 
any integer value in [0,ci]. Denote Di = {d, 0} as the weight 
pattern in which all d weights are not relayed. Then the weight 
pattern in general has one of two forms, Di = {d, 0 } and 
Dd 7 ^ Di. Using the MGF, the UPEP can be computed for 
general cases as follows: 


• Case 1: = Di. In this case, all d weights locate in 

the source-destination block, resulting in 72 = d'ysD and 
= 'I'tsd {ds). In this case we have: 


PuUlDi) 


1 

TT 

1 

2 


sin 

Jo sinP + d^sjj 



dlSD ^ 
1 -I- d'ygjj j 


(A.l) 


Case 2: 7 ^ Di. There is always d 2 weights are 

relayed, resulting in T'.y,,(s) = {ds) x vp-,,,, {d 2 s). 

Prom (0 we have: 

Nr Nr 

Pu(d|Drf)= ^((-1)^-1 ^ 


i=i 


1 


t/2 


sinf 


TT J (sin -h djso) (sin 0^ + d2lsei,j) 


de 


Nr Nr 

X] d:i{dTysr>,d2lSel,i)), 


i=i 


ni=l,...,nj=l 


(A.2) 


where 


Ii{aM = ^ M - “ 


a — b\ \ a h — ay 1 -\- b 


Appendix B 
Proof of Theorem|2] 

The diversity order is defined as the negative exponent of 
UPEP in log-scale when the average SNR 7 tends to infinity 


Td = - Jim 

7—>-oo 


/ logPu(d|Dd) \ 
V log7 J 


(B.l) 


Using the upper bound of UPEP p2\ as Pu (dlD^;) < 
|T'-y,,(l/2) < T'.y,,(l/2) and recall ( jB.!) ! we have 


Td> - Jim 

7—>-oo 


log ^71: (1/2) 

log 7 


(B.2) 
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Similar to Appendix we consider two cases. 

• Case 1; = Di. There is not any symbol helped by the 

relay and (1/2) = (<^/2). The diversity order 

in this case is given by 


Tfi > — lim 

7—J-OO 

> — Jim 

7—>-oo 


log id/2) 

log 7 

(1 + djsDl 2) 

log 7 


= 1 . 


(B.3) 


Then the UPEP has diversity order of 1 when ^2 = 0 and 
we can write Pu ((i|Di) A SNR“^. 

• Case 2; J Di. The MGF of 72 in this case has 
a form of (1/2) = (d/2) x (^ 2 / 2 ). 

Consequently, the diversity order is given as follows: 

Td> - lim — {d/2) log^Tsei {d2/2) 

7— >-oo log "7 7— >-oo log 7 

^ ■v' 

TSel 

= 1 + TSeh (B-4) 


where rsei is the diversity order of the best relay signal 
(without direct link). It has shown in ||^ that the best relay 
selection achieves diversity order is equal to the total 
number of available relays, we have rsei = Nr- There¬ 
fore the system diversity order in this case is equal to 
Nr +1. In order words, we can write Pu (d|Dd ^ Di) A 

Combine the two cases above we complete the proof of 
Theorem |2l 


Appendix C 
Proof of Theorem[3] 

It is noted that the total SNR in NCC is given by 7Sjvc = 
di^SiD + d 2 js 2 D + dfi^NC- To derive the UPEP for NCC, 
the MGF method is employed. Lemma [T] and Lemma state 
that there are at least two weights in {di,d 2 ,dii} are non¬ 
zero. Therefore, the weight pattern can only has one of 
these four cases: 1) da — 0, 2) di = 0, 3) d 2 — 0 and 4) 
did 2 d^ 0. Consequently, the total SNR ^y.mc one of 
corresponding four values. 

• Case 1: d^ = 0, there is not any weight on the relay 
channel. The total SNR has a form of ^y.nc ~ di'js-i_D + 
d 2 lS 2 D, and its MGF is given as follows: 

(dis) X 4 / 7320 (^ 25 ) 


^ - di'yS idS ^ - d2l s2ds' 

Now the UPEP Pu (dlW^) can be computed using the 
MGF method as follows: 

7r/2 

P. UIW,) = i / (^) 

0 

= Ti (rfi7siD> '^27 s2d) : (C.2) 

where Xi (a, b) has been dehned in Appendix [a| 

• Case 2: di = 0. In this case, the total Sl^ equals 

7Snc = d 2 'ys 2 D + dn'^NC- Given the PDF of 7^0 


in (doji, the MGF of the total SNR is given as follows: 

4/7S2£,(d2s) X (d^s) 

Nr Nr ^ ^ 

= yf(_i)/-i y -L-L-) 

""^275.05 1 -dR-fNC,jsJ 

(C.3) 

The UPEP Pu^<^|W(i^ is computed using the MGF 
method as follows: 

Nr Nr 


p,(diw,)=£((-iy-’ £ 


2 = 1 


ni<-"<ni 


t/2 


smf 


-dd 


equals 

Nr 


( sin 02 + d27s,B) ( sin + dnjjvcj) 

Nr- 

H((-l)^”^ ^i{d27s2D,dR7NC,j))- 

(C.4) 

,(d|Wd 


TV, 






• Case 3: do = 0. Similar to Case 2 we have P, 


Nr Nr 

^((-1)-^ ^ ^ Ti(^di7_5^£,,dij7jvc,j)y 


i=i 




Case 4: did 2 dR J 0. In this case, there are three weights 
in the total SNR, resulting in 7 ejvc = i^i 7 SiD+d 27 S 2 D + 
dR7NC- The MGF of 7s„c is given as follows: 

4/7s^£,(dis) X (d2s) X 4/-y^p(dijs) 


Nr Nr 

E(<-i)''‘ E 


i=i 


1 


1 


)■ 


1 - dllStOS 1 - d2ls2DS 1 - dRjNcjS 

(C.5) 

Appying the MGF method to compute the UPEP, we 
have: 


Nr Nr 

Pu(<i|wj)=y;((-iy-' y; 


J =1 




1 

TT 


Tr/2 


nir^---ytnj 


Sin 




-dd 


(sin d2 -f di 7 s^ £,) (sin 0^ + d27s2D) (sin d2 -f dR^j^c.j) 

Nr- Nr 

= d:2{di7siD^d27s2D,dR7ffcj 

j — 1 


(C. 6 ) 


where l 2 {a,b,c) has been dehned in Theorem]^ 
Combining four cases above gives Theorem 
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Appendix D 
Proof of Theorem|4] 

Similar to Appendix we employ the upper bound of 
UPEP to derive diversity order for NCC; 


T„ > — _lim 

-y—^oo 




log 7 


- ■ 


(D.l) 


We consider four cases: 

• Case 1: dij = 0. In this case, all weights locate in the 
source-destination channels, resulting in ip. 




(l/2) 

X ^(^ 2 / 2 ^. The diversity order in 


this case is given by: 

log 

log 7 

(1 + 


Tn > — _lim 

7—>-oo 


— Jim 

7—>-oo 

-1 


l0g^7S2D(c^2/2) 


log 7 


> - Jim 

7—>-oo 


— Jim 

7—>-oo 

= 1 + 1 = 2 . 


log 7 

(1 + d2lS2D/‘^ 


-1 


log 7 


(D.2) 


This is enough to say the UPEP has diversity order of 

= o'! = 


u{d\dR 


2 when dn = 0 and we can write P 
SNR- 2 . 

• Case 2: di = 0. The MGE of the total SNR in this 


case has a form of T*.- 


(l/2 


= Ti 


■yS2D 


(d 2 / 2 ) 


'I'tnc {ydfi/2^. Consequently, the diversity order is given 


as follows: 


> - Jim 

7—>-oo 


= 1 + J, 


logT' 


732 0 


(^ 2 / 2 ) 


log 7 


— Jim 

7—>-oo 


logT-^ 


■ {dR/2 


log 7 


J 


(D.3) 


where J is the diversity order of the best relayed signal 
without the direct link, which equals diversity order of 
the best relay selection for two-way relay channels. It 
has shown in ||^ that this diversity order is equal to Nr- 
Therefore we have the system diversity order when di = 

0 is equal to AJ + 1. In order words, Pu|^d|(ii = 0^ = 
SNR-^’-\ 

• Case 3: ^2 = 0. Similar to case 2 we have the diversity 
order equals AJ + 1 . 

• Case 4: did 2 dR ^ 0. In this case the MGE of 7 Ejvc 1^ 
a product of three terms: 

(C3) 

= »™.o(A/2) X J/ 2 ) X J«/2). 

Substituting ( |D.4[ ) into fDT) we have 

log^'7Sio('^i/2) logT'^s^o(d 2 / 2 ) 

Tn>- lim -^-+- - - lim -^—+- - 

7 —>-oo log 7 7 ->-oo log 7 


— lim 

7—>-oo 


log 

^7AfC (rffl/ 2 ) 


log 7 


J 

= 1 + 1 + = Nr T 2. 


(D.5) 


We can write P 


u^d|did2d_R 7^ 0^ A SNR 


-Nr-2 


Prom four cases above we have Theorem proved. 
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